Astronomy & Astrophysics manuscript no. 7605 


© ESO 2008 


February 1, 2008 





An ultraluminous supersoft source with a 4 hour modulation in 

NGC 4631 

S. Carpano\ A. M. T. Pollock\ A. R. King^, J. Wilms^ and M. Ehle' 



o 
o 



o 



> 

m 
m 

p 
O 

o 



1 XMM-Newton Science Operations Centre, ESAC, ESA, PO Box 50727, 28080 Madrid, Spain 
' Department of Physics and Astronomy, University of Leicester, Leicester LEI 7 RH, UK 

^ Dr. Remeis-Observatory, Astronomisclies Institut der FAU Erlangen-Niirnberg, Sternwartstr 7, 96049 Bamberg, Germany 



Submitted: 4 April 2007; Accepted: 6 July 2007 



ABSTRACT 



Context. Supersoft X-ray sources (SSSs) are characterised by very low temperatures (<100 eV). Classical SSSs have bolometric 
luminosities in the range of 10^*- 10^* ergs" 'and are modelled with steady nuclear burning of hydrogen on the surfaces of white 
dwarfs. However, several SSSs have been discovered with much higher luminosities. Their nature is still unclear. 
Aims. We report the discovery of a 4 h modulation for an ultraluminous SSS in the nearby edge-on spiral galaxy NGC 4631, observed 
with XMM-Newton in 2002 June. Temporal and spectral analysis of the source is performed. 

Methods. We use a Lomb-Scargle periodogram analysis for the period search and evaluate the confidence level using Monte-Carlo 
simulations. We measure the source temperature, flux and luminosity through spectral fitting. 

Results. A modulation of 4.2 + 0.4 h (3cr error) was found for the SSS with a confidence level >99%. Besides dips observed in the 
light curve, the flux decreased by a factor of 3 within ~ 10 h. The spectrum can be described with an absorbed blackbody model 
with kT ~ 67 eV. The absorbed luminosity in the 0.2-2 keV energy band was 2.7 x 10^** ergs"' while the bolometric luminosity was 
a hundred time higher (3.2 x lO""' erg s"'), making the source one of the most luminous of its class, assuming the best fit model is 
correct. 

Conclusions. This source is another very luminous SSS for which the standard white dwarf interpretation cannot be applied, unless 
a strong beaming factor is considered. A stellar-mass black hole accreting at a super Eddington rate is a more likely interpretation, 
where the excess of accreted matter is ejected through a strong optically-thick outflow. The 4 h modulation could either be an eclipse 
from the companion star or the consequence of a warped accretion disk. 
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1. Introduction 



The X-ray binaries known as supersoft sources (SSSs) are char- 
acterised by very soft emission with temperatures <100eV and 
bolometric luminosities exceeding lO^'^ergs"'. The standard 
model for such soft and h igh luminosities was proposed by 
. Ivan den Heuvel et~an d 19921) as nuclear burning of hydrogen on 
' the surface of a white dwarf of mass in the range 0.7-1.2Mo. 
' On the other hand, if luminosities greatly exceed the Eddington 
limit for solar mass compact objects, it is not obvious that 
this type of model could apply unless beaming is considered. 
Several such very lumin ous SSSs have been reported in nearby 



galaxies : two in M 101 (iDi Stefano & Kong]|2 003!: 'Ko ng et al 
2004; Mukai et al.1 120051). one in M 51 jPi Stefano & Kong 



2003), o ne in M 81 dSwartz et al.ll2002h . o ne in the Antennae 
(Fabbiano et al.l2003b and two in NGC 300 jKong & Pi Stefanol 
120031; ICarpano et al.ll2006l) . 

In one of these SSSs, a short period of 5.4h has been 
reported: the source located in the face-on spiral galaxy 
NG C 300, which was di scussed by iKong & Pi Stefanol j2003l) 
and lCarpano et al.l (|2006|) . The modulation was only present dur- 
ing a single XMM-Newton observation but was not visible 6 
days earlier, when the source was a few times brighter. In this 
Letter, we report the discovery of an even more luminous SSS 
with a slightly shorter modulation, located in the nearby edge- 
on galaxy NGC 4631. 
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NGC 4631 is a SBc/d t ype galaxy located at an as- 
sumed distance of 7.5 Mpc (iGolla & Hummell 11994 . with 
low galactic foreground absorption (A^h - 1-2 x 10^° cm"^; 
Pickey & Lockm an 1 990)) . A study of the giant di ffuse X-ray 
emittin g corona around the g alaxy was performed b v lWang et al.l 
(1200 lb with Chandra and by lTiillmann et alJ (l2006l) with XMM- 
Newton. The point- source popula tion viewed by ROSA T was 
studied bv lVogler & Pietschid 19961) and lRead et al](ll997l) . They 
reported the detection of 7 sources, including the supersoft 
source described in this Lett er. Spectral fitting of this SSS us- 
ing a bremsstrahlung model dRead et al.lll997h gave a tempera- 
ture of 0.1 keV with high (8 x 10^' cm"^) intrinsic absorption; 
the corresponding absorbed luminosity in the 0. 1-2.0 keV en- 
ergy band was of 1 .9 x 10^'^ erg s"'. In the long-term light curve 
(IVogler & Pietschlll996h the source was not visible in two of six 
data sets although, because the source flux was close to the de- 
tection limit in the other four, time variability could not be es- 
tablished. 



In this Letter, we report the discovery of a 4h modulation 
for the ultraluminous SSS in NGC 4631. The remainder of the 
Letter is organised as follows. Section |2] describes the XMM- 
Newton observation and data reduction. In Sect. [3] we present 
timing and spectral analysis of the SSS, reserving discussion for 
Sect.!] 
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Fig. 1. Mean XMM-Newton EPIC-MOS and pn background- 
subtracted light curve of the SSS in NGC 463 1 . Error bars are at 
1 cr. The straight line shows the best linear fit to the lightcurve. 

2. Observation and data reduction 

NGC 4631 was observed by XMM-Newton on 2002 June 8 
for 55ksec. T he EPIC-MOS jTurner et al.ll200lh and EPIC-pn 
dStriider et al.i 2001) cameras were operated in full frame mode 
with the medium MOS filter and the thin pn filter After screen- 
ing the MOS data for proton flares using standard procedure^ 
a total of 45 and 37 ksec of low-background emission remained 
for the MOS and pn, respectively. 

Using the XMM-Newton Software Analysis System (SAS) 
edetect_chain task, which performs maximum-likelihood 
source detection, the SSS was detected with a maximum 
likelihood of 1919. The best-fit coordinates were: ffjaooo - 
12H2™16^1 and ^j^ooo = +32°32'49'.'5 with a statistical error 
ofO'.'3. 



3. Timing and spectral analysis of the SSS 

The combined EPIC-MOS and pn background-subtracted 
lightcurve is shown in Fig. [1] with a time bin size of 500 sec. 
Times are given in hours from the start of the observation. 
Periods of high background at the end of the observation have 
been excluded from the data. Two dips are apparent in the 
lightcurve in addition to a long-term decrease by a factor of 3 
(Fig- [II straight line). 

We performed a period search with Lom b-Scargle peri- 
odogram analysis (lLomblll976l: IScarglelll982h after having re- 
moved the general trend of the decreasing flux. Confidence lev- 
els were measured by means of Monte Carlo simulations assum- 
ing a null hypothesis of white noise. Results are shown in Fig. [J] 
We found that the 4.17h period is significant at > 99% confi- 
dence. This modulation, although less significant, is also present 
in the MOS data alone. We estimated the uncertainty of the pe- 
riod by fitting a sine function using the IDlE task curvefit to 
the detrended light curve. The best-fit period by this method was 
4.11^[j!|^h, close to the value from the Lomb-Scargle analysis. 
The errors quoted are at the 3cr level. 

The corresponding light curve folded at 4.17 h is shown in 
Fig. [3] The shape is slightly asymmetric, the flux change before 
the dip sharper than afterwards. The dip is also relatively narrow. 




Period (hours) 

Fig. 2. Search for periodicities in the XMM-Newton light curve 
of the NGC 4631 SSS using Lomb-Scargle periodogram analy- 
sis, after removing the decreasing flux trend. The full, dotted and 
dashed lines represent the 68%, 90% and 99% confidence levels, 
respectively, as determined from Monte Carlo simulations. 
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Fig. 3. XMM-Newton light curve of the NGC 4631 SSS folded 
at 4. 17 h using 30 bins, after removing the decreasing flux trend. 

Fig.|4]shows the pn and MOS spectra of the SSS. The data 
are binned to have at least 25 counts in each energy bin. Almost 
no flux is visible above 2keV. We first tried to model the data 
with an absorbed single blackbody, yielding = 1.95. From 
the residuals shown in Fig. |4] some excess around 0.85 keV is 
observed. Adding a Gaussian line significantly improved the fit, 
resulting in a reduced - 1.22 for v = 53 degrees of freedom . 
The absorbed bremsstrahlung model used bv lRead et al.ld 19971) . 
with or without a Gaussian line, does not describe the data well 
(xl - 3). We also tried to fit a black body with a power-law to 
describe the continuum but the normalisation constant of this 
last component is very small and spectral parameters as well 
as luminosities do not change significantly (less than a factor 
of 2). Alternatively to the gaussian emission line at 0.8 keV, an 
absorption edge could also describe the spectrum (xl - 1.35). 
Unfortunately due to the low statistic of the data above 1 keV, the 
' http : //ximti. esac . esa . int/ external/ximn_user_support/docu mbscpittlon/daptliii^flaSfi^ler cannot be constrain ed. Absorption 
^ [http : //www ■ ittvis ■ com/idl| edges have been reported bv lSwartz et al.l (|2002|) and lKong et al.l 
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Table 1. Results of the spectral fits for the SSS in NGC 4361, 
using an absorbed blackbody model and a Gaussian line 
(phabs* (bbody+gauss), in XSPEC). 



Parameter Value (units) 



A'h 


4.94![ ^^ X 10-'(cm-^) 




kT 


61% (eV) 




El 


0.80+°;°3 (keV) 






0.07+«;«* (keV) 




-fo.2-2keV 


3.97+f J* X 10-'-* (ergs 




I obs 

'"0.2-2keV 


2.67!3-20xl038(ergs- 


) 


J unabs 
'"0.2-2keV 


1.99+']^^ X 10'"' (ergs- 


) 


^'bol 


3.22+^35 X 10"° (ergs- 


) 



Fig. 4. XMM-Newton EPIC-pn and MOS spectra of the SSS in 
NGC 4361 and the best-fit spectral model with the apparent line 
emission. Bottom: residuals expressed in cr. 



(|2004 for the luminous SSSs in M 81 and M 101, respectively, 
indicating the presence of a warm absorber in the vicinity of the 
X-ray source. The best-fit parameters of the absorbed combined 
blackbody and single Gaussian-line model are shown in Table[T] 
A^H is the equivalent column density of neutral hydrogen, kT the 
temperature, E]^ the energy of the line and ctl its width. The 
corresponding 0.2-2 keV flux, absorbed, unabsorbed and bolo- 
metric luminosities are shown in the last four rows. Uncertainties 
are given at the 90% confidence level, except for the bolometric 
luminosity at 99%. As the source was close to a pn CCD gap, the 
flux and luminosities were calculated from the MOS data alone. 

Our data modelling shows the source to have a very high 
intrinsic absorption column giving a bolometric luminosity that 
makes the SSS one of the most luminous of its class. Such high 
absorption could be due to a strong intrinsic absorption in the 
vicinity of the emitting source or due to gas of the hosting galaxy 
since it is edge-on oriented. The source is also highly variable. 
The object has also been observed once with Chandra, about 
two years before the XMM-Newton observation, on 2000 April 
16 for 59ksec when 31 net counts were detected. Converting 
this rate to flux using WebPIMM^ with the spectral parameters of 
Table [T] we estimate that the mean observed luminosity at that 
time in the 0.2-2 keV energy band was 1 .7 x lO-'^ erg s"' . This is 
about a factor of 10 lower than the XMM-Ne wton measuremen t 
two years later and the earlier ROSAT value (iRead et al.lll997h . 

4. Discussion 

The SSS reported in this Letter is one of the most luminous of 
its class. Besides the bright and soft X-ray emission, the source 
is particular for its 4 h modulation superimposed on a decreas- 
ing flux trend during the XMM-Newton o bservation. On only 
one o ther occasion, for a SSS in NGC 300 jKong & Di Stefanol 
l2003h . has a short modulation been reported for such a bright 
SSS, of 5.4 h in that case, in the one of the two XMM-Newton 
observations in which the source was fainter A discuss ion of the 
source has also been reported bv lCarpano et al.l ( |2006|) . 

The high luminosity (> 10^" erg s"') of the NGC 4361 SSS, 
which exceeds the Eddington luminosity for a white dwarf or 
a neutron star by a factor of 100, might be due to a strong 



wind or outflow of a stellar-mass black hole accreting at a 
super Eddington rate . Th is idea has also b een suggested by 
iFabbiano etal] (|2003|) and lMukai eTal] (l2005b for flie ultralumi- 
nous SSSs in the Antennae and M 101, respectively. The mode l 
has been reconsidered in more detail bv lKing & Poundsl ( l2003h . 
as summarised below. Other less convincing models in the litera- 
ture, as summarised bv lCarpano et al.l (l2006h . have variously in- 
volved strong beaming from steady nuclear burning of hydrogen 
at the surface of a white dwarf; accretion onto an intermediate- 
mass black hole; or soft emission of an accretion disk observed 

at high inclination. 

lKing& Poundsl (120031) suggested that a black hole accreting 
at super-Eddington rate, ejects the excess matter through a strong 
wind or outflow that is likely to be optically thick to electron 
scattering for Mom ~ ^Edd with a photospheric radius of the 
order of few tens of the Schwarzschild radius of the black hole. 
The outflow is radial and occupies a double cone of solid angle 
Anb, where the beaming factor h can be either ~ 1 or «: 1 . The 
emission line present in the spectrum of the SSS in NGC 4631, 
which is likely to be composed of several unresolved lines, or 
alternatively the absorption edge, might come from optically thin 
parts of the matter outflow. These could either be near the edges 
of the conical outflow, or the leading edge of the outflow if this is 
time-depen dent. The temperature resulting from the black body 
enussion is (iKing & Poundsll2003h : 



(1) 



where g(b) = 1/^ or l/(2b^^^) (for ~ 1 or ^ « 1), Mi = 
Mout/(lM0yr-'), and Ms = M/IO^^Mq, and M is flie mass of 
the accretor 

Assuming the source luminosity is close to the Eddington 
limit, Mout ~ M - L/irjc^), where L is the luminosity and t] the 
radiative efficiency, we are able to estimate the mass of the black 
hole 



M = 0.27 



,1/4 



4/3 



Mr 



(2) 



|http : //heasarc . gsf c . nasa . gov/Tools/wBpimms . html | 



As g''''* is of order unity unless b < 0.1, for a typical value of 
T] ~ 0.1, the bolometric luminosity of 3.2 x 10'*" erg s"' gives a 
corresponding black-hole mass of 6M0 and an accretion rate 
of 5.6 X 10"^ Moyr-'. Such an accretion rate can easily be 
achieved either by thermal-time-scale mass transfer after a high- 
mass companion has filled its Roche Lobe or during outbursts of 
soft X-ray transients (King 2002). 

The 4h modulation could be caused either by an eclipse 
from the companion star or by the presence of a warped disk. 
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This latter hypothesis is driven by the fact that the 5.4 h modu- 
lation in the SSS of NGC 300 was observed in only one of two 
XMM-Newton observations. Barnard et al. (200^ have recently 
reported that the X-ray binary Bo 158 in M 31, shows regular 
dips at 2.78 h but not in all of their observations. These authors 
suggested this dipping to be caused by an elongated precessing 
disk. 

To summarise, the ultraluminous SSS in NGC 4631 reported 
in this Letter may be an X-ray binary system with a short (<10h) 
orbital period, where the high mass transfer rate is either due to a 
massive companion filling its Roche Lobe or because the system 
is in outburst. When matter is accreted above the Eddington rate, 
the excess matter is ejected via a strong, optically-thick outflow 
in which case a mass of 6 Mq can be estimated for the black-hole 
companion. 

Acknowledgements. This paper is based on observations obtained with XMM- 
Newton, an ESA science mission with instruments and contributions directly 
funded by ESA Member States and NASA. 



References 

Barnard, R., Foulkes, S. B., Haswell, C. A., Kolb, U., Osborne, J. P. & Murray, 

J. R. 2006, MNRAS, 366, 287 
Carpano. S., Wilms. J., Schirmer, M., & Kendziorra, E. 2006, A&A, 458, 747 
Dickey, J. M. & Lockman, R J. 1990, ARA&A, 28, 215 
Di Stefano, R. & Kong, A. K. H. 2003, ApJ, 592, 884 

Fabbiano, G., King, A. R., Zezas, A., Ponman, T. J., Rots, A. & Schweizer, F. 

2003, ApJ, 591,843 
Golla, G. & Hummel, E. 1994, A&A, 284, 777 
King, A. R. 2002 MNRAS, 355, L13 
King, A. R. & Pounds, K. A. 2003, MNRAS, 345, 657 
Kong, A. K. H. & Di Stefano, R. 2003, ApJ, 590, L13 
Kong, A. K. H., Di Stefano, R. & Yuan, R 2004, ApJ, 617, L49 
Lomb, N. R, 1976, Ap&SS, 39, 447 

Mukai, K., Still, M., Corbet, R. H. D., Kuntz, K. D. & Barnard, R. 2005, ApJ, 
634, 1085 

Read, A. M., Ponman, T. J. & Strickland, D. K. 1997. MNRAS, 286, 626 
Scargle, J. D., 1982, 263, 835 

Striider, L., Briel, U., Dennerl, K., et al. 2001, A&A, 365, L18 
Swartz, D. A., Ghosh, K. K., Suleimanov, V., Tennant, A. F. & Wu, K. 2002, 
ApJ, 574, 382 

TuUmann, R., Breitschwerdt, D., Rossa, J., Pietsch, W. & Dettmar, R.-J. 2006, 
A&A, 457, 779 

Turner, M. J. L., Abbey, A.. Arnaud, M., et al. 2001. A&A, 365, L27 

van den Heuvel, E. P. J., Bhattachai'ya, D., Nomoto, K. & Rappaport, S. A. 1992, 

A&A, 262, 97 
Vogler, A. & Pietsch, W. 1996, A&A, 311, 35 

Wang, Q. D., Immler, S., Walterbos, R., Lauroesch, J. T. & Breitschwerdt, D. 
2001, ApJ, 555, L99 



List of Objects 



'M 101' on pagejl] 
'M 51' on page 1 
'MSr onpageG 




